3d wireless optical positioning method and system

ABSTRACT

The present invention provides a 3D wireless optical positioning method and system, including the steps of: arranging two LED lamps on the ceiling to transmit optical information and provide illumination; arranging a receiver including two photodetectors in a receiving plane; calculating the distance between the LED lamps and the photodetectors respectively through the TOA (Time of Arrival) method; and finally determining the actual position and orientation angle of the receiver based on the geometrical relationship between the LED lamps and the photodetectors in the XYZ coordinate system, the two photodetectors having a distance determined as l therebetween and being situated in the same receiving plane, the receiver being situated below the two LED lamps, the range where the receiver is to be positioned being on any side of the plane consisting of the two LED lamps and the origin.

This application is the National Stage Application of PCT/CN2021/116639, filed on Sep. 6, 2021, which claims priority to Patent Application No. 202111016290.5, filed on Aug. 31, 2021, which is incorporated by reference for all purposes as if fully set forth herein.

FIELD OF THE INVENTION

The present invention is related to the technical field of positioning, and more particularly to a 3D wireless optical positioning method and system.

DESCRIPTION OF THE RELATED ART

Positioning demand is present everywhere in life. Traditional outdoor positioning techniques, such as Beidou and GPS, have meter-level positioning accuracy outdoors, but their signals will be severely attenuated when they pass through walls and other obstructions, so Beidou and GPS are not suitable for indoor positioning that requires high accuracy. In recent years, as an alternative to GPS, many new indoor positioning techniques, such as ultrasonic, Bluetooth, WiFi, visible light positioning and so on, have been continuously developed. Because of its advantages of energy saving, high accuracy, fast positioning speed, strong anti-interference and low cost, visible light positioning stands out among many indoor positioning techniques. However, due to the limitation of its positioning principle, there is still great room for improvement in 3D positioning and orientation, which limits the large-scale use of visible light positioning to a certain extent.

The common visible light positioning solutions have the following problems: in use of the scenario analysis method, the characteristic parameters of the environment where the LED light source is located need to be measured first. If the space of the area to be measured is large, the parameter measurement process requires a great amount of work. Secondly, when the environment changes, the parameters need to be re-determined, and the portability of scenario analysis method is poor.

In the traditional geometric measurement method, at least three LED light sources are needed in order to realize trilateration or triangulation, which limits the use of positioning solutions in scenarios with insufficient number of LEDs. In addition, with this method, it is difficult to determine the orientation of the receiving terminal because there is only one photoelectric detector at the receiving terminal.

In the image sensor method, because it deals with the change in projection of the light source in the image sensor, the 3D spatial information becomes the 2D plane information, and the image sensor method cannot directly calculate the height of the receiving terminal, but can only realize 2D positioning.

Therefore, current indoor visible light positioning methods are mainly intended for indoor 2D positioning, which makes it difficult to achieve accurate indoor 3D positioning, and also makes it difficult to determine the actual orientation of the receiving terminal.

SUMMARY OF THE INVENTION

It is an object of the present invention to provide a 3D wireless optical positioning method and system that enables 3D positioning and orientation of a terminal while providing illumination and enables accurate indoor 3D positioning and accurate determination of the actual orientation of the terminal even in the case of a reduced number of optical transmitters.

To address the technical problem mentioned above, the present invention provides a 3D wireless optical positioning method including the steps of:

-   -   arranging LED lamps on the ceiling to transmit optical         information and providing illumination and arranging a receiver         on a receiving plane to receive the optical information, in         which the LED lamps include a first LED lamp and a second LED         lamp with a coordinate of (x_(t1), y_(t1), z_(t1)) and (x_(t2),         y_(t2), z_(t2)) respectively, a first photodetector and a second         photodetector with a coordinate of ({circumflex over (x)}_(r1),         ŷ_(r1), {circumflex over (z)}_(r1)) and ({circumflex over         (x)}_(r2), ŷ_(r2), {circumflex over (z)}_(r2)) respectively are         arranged on the receiver, the distance between the first         photodetector and the second photodetector is defined as l,         meanwhile the first photodetector and the second photodetector         both face upwards on the receiving plane, the middle point         between the first photodetector and the second photodetector         defines the actual position of the receiving terminal to be         predicted, and the direction from the first photodetector to the         second photodetector defines the orientation of the receiving         terminal, that is, the included angle between the line         interconnecting the first photodetector and the second         photodetector and the positive half of the X axis is the         orientation angle η;     -   through the TOA (Time of Arrival) principle, measuring the time         required for the optical signal to be transmitted from the first         LED lamp and the second LED lamp to and received by the first         photodetector and the second photodetector respectively, and         multiplying the propagation time by the speed of light to         calculate the distance d₁ between the first LED lamp and the         first photodetector, the distance d₂ between the second LED lamp         and the first photodetector, the distance d₃ between the first         LED lamp and the second photodetector and the distance d₄         between the second LED lamp and the second photodetector; and     -   obtaining the actual position and orientation angle of the         receiver based on the geometrical relationship between the LED         lamps and the receiver in the XYZ coordinate system, the first         photodetector and the second photodetector having a distance         determined as l therebetween and being situated in the same         receiving plane, the receiver being situated below the first LED         lamp and the second LED lamp, the range where the receiver is to         be positioned being on any side of the plane consisting of the         first LED lamp, the second LED lamp and the origin.

As a further improvement of the present invention, the corresponding photodetector and LED lamp are in time synchronization in calculating the distances d₁-d₄.

As a further improvement of the present invention, the following equation set can be obtained for the distances d₁-d₄ in the XYZ coordinate system:

d ₁ ²=({circumflex over (x)} _(r1) −x _(t1))²+(ŷ _(r1) −y _(t1))²+({circumflex over (z)} _(r1) −z _(t1))²  (1)

d ₂ ²=({circumflex over (x)} _(r1) −x _(t2))²+(ŷ _(r1) −y _(t2))²+({circumflex over (z)} _(r1) −z _(t2))²  (2)

d ₃ ²=({circumflex over (x)} _(r2) −x _(t1))²+(ŷ _(r2) −y _(t1))²+({circumflex over (z)} _(r2) −z _(t1))²  (3)

d ₄ ²=({circumflex over (x)} _(r2) −x _(t2))²+(ŷ _(r2) −y _(t2))²+({circumflex over (z)} _(r2) −z _(t2))²  (4)

Meanwhile, the following supplementary equations can be obtained as the first photodetector and the second photodetector have a distance determined as l therebetween and are situated in the same receiving plane:

l ²=({circumflex over (x)} _(r2) −{circumflex over (x)} _(r1))²+(ŷ _(r2) −ŷ _(r1))²+({circumflex over (z)} _(r2) −{circumflex over (z)} _(r1))²  (5)

{circumflex over (z)} _(r2) ={circumflex over (z)} _(r1)  (6)

Given the known d₁, d₂, d₃ and d₄ and the formulas (5) and (6) and in combination with the fact that the receiver is situated below the first LED lamp and the second LED lamp and the range where the receiver is to be positioned is on any side of the plane consisting of the first LED lamp, the second LED lamp and the origin, the actual position of the receiver (({circumflex over (x)}_(r1)+{circumflex over (x)}_(r2))/2, (ŷ_(r1)+ŷ_(r2))/2, ({circumflex over (z)}_(r1)+{circumflex over (z)}_(r2))/2) and the orientation angle of the receiver can be obtained through solution of:

$\begin{matrix} {\hat{\eta} = {{{sign}\left( {{\hat{y}}_{r2} - {\hat{y}}_{r1}} \right)}{arc}{{\cos\left\lbrack \frac{\left( {{\hat{x}}_{r2} - {\hat{x}}_{r1}} \right)}{\sqrt{\left( {{\hat{x}}_{r1} - {\hat{x}}_{r2}} \right)^{2} + \left( {{\hat{y}}_{r1} - {\hat{y}}_{r2}} \right)^{2}}} \right\rbrack}.}}} & (13) \end{matrix}$

As a further improvement of the present invention, the process of solving the actual position of the receiver (({circumflex over (x)}_(r1)+{circumflex over (x)}_(r2))/2, (ŷ_(r1)+ŷ_(r2))/2, ({circumflex over (z)}_(r1)+{circumflex over (z)}_(r2))/2) and the orientation angle of the receiver specifically includes the following steps:

-   -   S1: as the formulas (1) and (2) each represent on the physical         sense a sphere, subtracting the formula (2) from the formula (1)         to obtain the equation for the plane P1 where the two spheres         intersect each other, the first photodetector being positioned         on the circle where this plane intersects the sphere represented         by the formula (1), this circle having the centre at K₁, the         coordinate of which being (a₁, b₁, c₁), and the radius of R₁;         and likewise, subtracting the formula (4) from the formula (3)         to obtain the equation for the plane P2 where the two spheres         intersect each other, the second photodetector being positioned         on the circle where this plane intersects the sphere represented         by the formula (3), this circle having the centre at K₂, the         coordinate of which being (a₂, b₂, c₂), and the radius of R₂,         the coordinate of the circle centre and the radius of the circle         being represented by the formulas:

K _(i)(a _(i) ,b _(i) ,c _(i))=(x _(t1)+(x _(t2) −x _(t1))w _(i) /L, y _(t1)+(y _(t2) −y _(t1))w _(i) /L, z _(t1)+(z _(t2) −z _(t1))w _(i) /L), i=1,2 . . .   (7)

R ₁=√{square root over (d ₁ ² −w ₁ ²)}

R ₂=√{square root over (d ₃ ² −w ₂ ²)}   (8)

where w_(i) represents the distance between the first LED lamp and the plane P_(i) (i=1, 2) where the two spheres intersect each other and L represents the distance between the first LED and the second LED lamp;

-   -   S2: performing two coordinate system transformations, where in         the first coordinate transformation, an XYZ coordinate system is         transformed into an X′Y′Z′ coordinate system, the X′Y′Z′         coordinate system having the circle centre K₁ as its origin of         coordinate, the straight line K₁ K₂ as the new Z′ axis, the         intersecting line of the circle K₁ and the original XY plane as         the new X′ axis; and in the second coordinate transformation,         the new X′Y′Z′ coordinate system is transformed into a         cylindrical coordinate system, whereupon after the two         coordinate transformations, the coordinates of the first         photodetector and the second photodetector are changed into (R₁,         Φ₁, 0) and (R₂, Φ₂, S) respectively, and meanwhile, in the         cylindrical coordinate system, the formula (5) and the         formula (6) are represented respectively as:

Φ₁−Φ₂=±arccos M  (9)

β(R ₂ sin Φ₂ −R ₁ sin Φ₁)=−γS  (10)

where M=(R₁ ²+R₂ ²+S²−l²)/(2R₁R₂), the distance between the two circle centres K₁ and K₂ is expressed as S=√{square root over (a²+b²+c²)}, β=−√{square root over ((a²+b²))}/S, γ=c/S, a=a₂−a₁, b=b₂−b₁, c=c₂−c₁;

-   -   S3: calculating Φ₁ and Φ₂ according to the formula (9) and the         formula (10), then performing two inverse coordinate         transformations to recover the coordinates of the first         photodetector and the second photodetector in the XYZ coordinate         system through the formulas (11) and (12), i.e.:

({circumflex over (x)} _(r1) ,ŷ _(r1) ,{circumflex over (z)} _(r1))′=(e _(x) ,e _(y) ,e _(z))(R ₁ cos Φ₁ ,R ₁ sin Φ₁,0)′+(a ₁ ,b ₁ ,c ₁)′  (11)

({circumflex over (x)} _(r2) ,ŷ _(r2) ,{circumflex over (z)} _(r2))′=(e _(x) ,e _(y) ,e _(z))(R ₂ cos Φ₂ ,R ₂ sin Φ₂ ,S)′+(a ₂ ,b ₂ ,c ₂)′  (12)

where

e_(x)=(b/√{square root over (a²+b²)}, −a/√{square root over (a²+b²)}, 0)′, e_(y)=[ac/(S√{square root over ((a²+b²))}), bc/(S√{square root over ((a²+b²))}), −√{square root over ((a²+b²)}/S]′, e_(z)=(a/S, b/S, c/S)′ are the orthogonal basis of the first coordinate transformation.

As a further improvement of the present invention, in the step S3, four sets of solution are obtained by solving Φ₁ and Φ₂, and accordingly, four sets of coordinates of the first photodetector and the second photodetector in the XYZ coordinate system are obtained through two inverse coordinate transformations, whereas the actual position includes only one set, and as the four sets of solution are spatially symmetrical with respect to the line interconnecting the first LED lamp and the second LED lamp, the real solution can be obtained through determination based on the following conditions, including specifically the following steps:

-   -   S4: excluding two sets of solution representing a position above         the first LED lamp and the second LED lamp considering the fact         that the position of the receiver is below the first LED lamp         and the second LED lamp;     -   S5: obtaining the single real solution out of the remaining two         sets of solution by restricting the receiver in movement on any         side of the plane consisting of the first LED lamp, the second         LED lamp and the origin;     -   S6: obtaining the coordinate of the receiver as (({circumflex         over (x)}_(r1)+{circumflex over (x)}_(r2))/2, (ŷ_(r1)+ŷ_(r2))/2,         ({circumflex over (z)}_(r1)+{circumflex over (z)}_(r2))/2)         through the single set of solution, in which the orientation         angle η of the receiver has a value expressed by the         formula (13) as:

$\begin{matrix} {\hat{\eta} = {{{sign}\left( {{\hat{y}}_{r2} - {\hat{y}}_{r1}} \right)}{arc}{\cos\left\lbrack \frac{\left( {{\hat{x}}_{r2} - {\hat{x}}_{r2}} \right)}{\sqrt{\left( {{\hat{x}}_{r1} - {\hat{x}}_{r2}} \right)^{2} + \left( {{\hat{y}}_{r1} - {\hat{y}}_{r2}} \right)^{2}}} \right\rbrack}}} & (13) \end{matrix}$

where sign is the sign function.

As a further improvement of the present invention, the number of the LED lamps is defined depending on the region where they are to be positioned.

A 3D wireless optical positioning systems includes: LED lamps including a first LED lamp and a second LED lamp with a coordinate of (x_(t1), y_(t1), z_(t1)) and (x_(t2), y_(r2), z_(t2)) respectively, arranged on the ceiling to transmit optical information and provide illumination;

-   -   a receiver provided with a first photodetector and a second         photodetector with a coordinate of ({circumflex over (x)}_(r1),         ŷ_(r1), {circumflex over (z)}_(r1)) and ({circumflex over         (x)}_(r2), ŷ_(r2), {circumflex over (z)}_(r2)) respectively,         arranged on a receiving plane to receive the optical         information;     -   in which the distance between the first photodetector and the         second photodetector is defined as l, meanwhile the first         photodetector and the second photodetector both face upwards on         the receiving plane, the middle point between the first         photodetector and the second photodetector defines the actual         position of the receiving terminal to be predicted, and the         direction from the first photodetector to the second         photodetector defines the orientation of the receiving terminal,         that is, the included angle between the line interconnecting the         first photodetector and the second photodetector and the         positive half of the X axis is the orientation angle η;     -   through the TOA (Time of Arrival) principle, the time required         for the optical signal to be transmitted from the first LED lamp         and the second LED lamp to and received by the first         photodetector and the second photodetector respectively is         measured, and the propagation time is multiplied by the speed of         light to calculate the distance d₁ between the first LED lamp         and the first photodetector, the distance d₂ between the second         LED lamp and the first photodetector, the distance d₃ between         the first LED lamp and the second photodetector and the distance         d₄ between the second LED lamp and the second photodetector; and     -   the actual position and orientation angle of the receiver are         obtained based on the geometrical relationship between the LED         lamps and the receiver in the XYZ coordinate system, the first         photodetector and the second photodetector having a distance         determined as l therebetween and being situated in the same         receiving plane, the receiver being situated below the first LED         lamp and the second LED lamp, the range where the receiver is to         be positioned being on any side of the plane consisting of the         first LED lamp, the second LED lamp and the origin.

As a further improvement of the present invention, the first LED lamp and the second LED lamp and the first photodetector and the second photodetector have synchronized operation time.

As a further improvement of the present invention, the following set of equations for the distances d₁ to d₄ in the XYZ coordinate system is obtained based on the geometrical relationship between the LED lamps and the receiver in the XYZ coordinate system:

d ₁ ²=({circumflex over (x)} _(r1) −x _(t1))²+(ŷ _(r1) −y _(t1))²+({circumflex over (z)} _(r1) −z _(t1))²  (1)

d ₂ ²=({circumflex over (x)} _(r1) −x _(t2))²+(ŷ _(r1) −y _(t2))²+({circumflex over (z)} _(r1) −z _(t2))²  (2)

d ₃ ²=({circumflex over (x)} _(r2) −x _(t1))²+(ŷ _(r2) −y _(t1))²+({circumflex over (z)} _(r2) −z _(t1))²  (3)

d ₄ ²=({circumflex over (x)} _(r2) −x _(t2))²+(ŷ _(r2) −y _(t2))²+({circumflex over (z)} _(r2) −z _(t2))²  (4)

Meanwhile, the following supplementary equations can be obtained as the first photodetector and the second photodetector have a distance determined as l therebetween and are situated in the same receiving plane:

l ²=({circumflex over (x)} _(r2) −{circumflex over (x)} _(r1))²+(ŷ _(r2) −ŷ _(r1))²+({circumflex over (z)} _(r2) −{circumflex over (z)} _(r1))²  (5)

{circumflex over (z)} _(r2) ={circumflex over (z)} _(r1)  (6)

Given the known d₁, d₂, d₃ and d₄ and the formulas (5) and (6) and in combination with the fact that the receiver is situated below the first LED lamp and the second LED lamp and the range where the receiver is to be positioned is on any side of the plane consisting of the first LED lamp, the second LED lamp and the origin, the actual position of the receiver (({circumflex over (x)}_(r1)+{circumflex over (x)}_(r2))/2, (ŷ_(r1)+ŷ_(r2))/2, ({circumflex over (z)}_(r1)+{circumflex over (z)}_(r2))/2) and the orientation angle of the receiver can be obtained through solution of:

$\begin{matrix} {\hat{\eta} = {{{sign}\left( {{\hat{y}}_{r2} - {\hat{y}}_{r1}} \right)}{arc}{{\cos\left\lbrack \frac{\left( {{\hat{x}}_{r2} - {\hat{x}}_{r1}} \right)}{\sqrt{\left( {{\hat{x}}_{r1} - {\hat{x}}_{r2}} \right)^{2} + \left( {{\hat{y}}_{r1} - {\hat{y}}_{r2}} \right)^{2}}} \right\rbrack}.}}} & (13) \end{matrix}$

As a further improvement of the present invention, the first LED lamp, the second LED lamp, the first photodetector and the second photodetector are provided with a time synchronization device.

The present invention has the following beneficial effects. In the present invention, a small number of optical transmitters is utilized and LED lamps are used as the light source, so that the deployment is simpler, the portability is better and the limitation is reduced. The present positioning method can be applied to various indoor scenarios to achieve accurate indoor 3D positioning without the need for additional devices, such as image sensors, and by use of the pair of photodetectors, the actual orientation of the terminal can be accurately determined. For the scenario of a row of lamps indoors, the method of the present invention is particularly applicable by positioning on one side of the lamps. That is, the present method has good extensibility and good portability, and can be stably used in various indoor scenarios.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic structural view of a system according to the present invention;

FIG. 2 is a schematic view showing intersection of two spheres in the step S1 according to the present invention;

FIG. 3 is a schematic view showing coordinate transformations in the step S2 according to the present invention;

FIG. 4 is a schematic view of 3D positioning results according to an embodiment of the present invention; and

FIG. 5 is the diagram of an accumulated distribution function of the estimation error for the orientation angle according to an embodiment of the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

The invention will be further explained with reference to the following drawings and particular embodiments, so that those skilled in the art can better understand and implement the present invention. However, the listed embodiments should not be taken as limitation of the present invention.

Referring to FIGS. 1-3 , an embodiment of the present invention provides a 3D wireless optical positioning method, including the following steps:

-   -   arranging LED lamps on the ceiling to transmit optical         information and providing illumination and arranging a receiver         on a receiving plane to receive the optical information, in         which the LED lamps include a first LED lamp and a second LED         lamp with a coordinate of (x_(t1), y_(t1), z_(t1)) and (x_(t2),         y_(t2), z_(t2)) respectively, a first photodetector and a second         photodetector with a coordinate of ({circumflex over (x)}_(r1),         ŷ_(r1), {circumflex over (z)}_(t1)) and ({circumflex over         (x)}_(r2), ŷ_(r2), {circumflex over (z)}_(r2)) respectively are         arranged on the receiver, the distance between the first         photodetector and the second photodetector is defined as l,         meanwhile the first photodetector and the second photodetector         both face upwards on the receiving plane, the middle point         between the first photodetector and the second photodetector         defines the actual position of the receiving terminal to be         predicted, and the direction from the first photodetector to the         second photodetector defines the orientation of the receiving         terminal, that is, the included angle between the line         interconnecting the first photodetector and the second         photodetector and the positive half of the X axis is the         orientation angle η;     -   through the TOA (Time of Arrival) principle, measuring the time         required for the optical signal to be transmitted from the first         LED lamp and the second LED lamp to and received by the first         photodetector and the second photodetector respectively, and         multiplying the propagation time by the speed of light to         calculate the distance d₁ between the first LED lamp and the         first photodetector, the distance d₂ between the second LED lamp         and the first photodetector, the distance d₃ between the first         LED lamp and the second photodetector and the distance d₄         between the second LED lamp and the second photodetector         respectively; and     -   obtaining the actual position and orientation angle of the         receiver based on the geometrical relationship between the LED         lamps and the receiver in the XYZ coordinate system, the first         photodetector and the second photodetector having a distance         determined as l therebetween and being situated in the same         receiving plane, the receiver being situated below the first LED         lamp and the second LED lamp, the range where the receiver is to         be positioned being on any side of the plane consisting of the         first LED lamp, the second LED lamp and the origin.

Specifically, it is proposed that the system consists of two elements as a transmitter and a receiver. The transmitter includes two LED lamps installed on the ceiling, LED1 and LED2, with a coordinate of (x_(t1), y_(t1), z_(t1)) and (x_(t2), y_(t2), z_(t2)) respectively, that can provide illumination. The plane where the receiver is situated is the receiving plane. A pair of (two) photodetectors is installed on the receiver. The two photodetectors are PD1 and PD2 respectively, with a coordinate of ({circumflex over (x)}_(r1), ŷ_(r1), {circumflex over (z)}_(r1)) and ({circumflex over (x)}_(r2), ŷ_(r2), {circumflex over (z)}_(r2)) respectively, configured to receive the optical information transmitted from the light source. The distance between PD1 and PD2 is l, and PD1 and PD2 both face upwards. The middle point between PD1 and PD2 defines the actual position of the receiving terminal to be predicted. The direction from PD1 to PD2 defines the orientation of the receiving terminal. The included angle between the line interconnecting PD1 and PD2 and the positive half of the X axis is defined as the orientation angle η.

Positioning: through the TOA (Time of Arrival) principle, the time required for the optical signal to be transmitted from LED1 and LED2 to and received by PD1 and PD2 respectively is measured, and the propagation time is multiplied by the speed of light to calculate the distance d₁ between LED1 and PD1, the distance d₂ between LED2 and PD1, the distance d₃ between LED1 and PD2 and the distance d₄ between LED2 and PD2. The following set of equations is obtained:

d ₁ ²=({circumflex over (x)} _(r1) −x _(t1))²+(ŷ _(r1) −y _(t1))²+({circumflex over (z)} _(r1) −z _(t1))²  (1)

d ₂ ²=({circumflex over (x)} _(r1) −x _(t2))²+(ŷ _(r1) −y _(t2))²+({circumflex over (z)} _(r1) −z _(t2))²  (2)

d ₃ ²=({circumflex over (x)} _(r2) −x _(t1))²+(ŷ _(r2) −y _(t1))²+({circumflex over (z)} _(r2) −z _(t1))²  (3)

d ₄ ²=({circumflex over (x)} _(r2) −x _(t2))²+(ŷ _(r2) −y _(t2))²+({circumflex over (z)} _(r2) −z _(t2))²  (4)

Meanwhile, the following supplementary equations can be obtained as PD1 and PD2 have a distance determined as l therebetween and are situated in the same receiving plane:

l ²=({circumflex over (x)} _(r2) −{circumflex over (x)} _(r1))²+(ŷ _(r2) −ŷ _(r1))²+({circumflex over (z)} _(r2) −{circumflex over (z)} _(r1))²  (5)

{circumflex over (z)} _(r2) ={circumflex over (z)} _(r1)  (6)

The equation set including formulas (1)-(6) is solved by the following steps:

First step: as shown in FIG. 2 , as the formulas (1) and (2) each represent on the physical sense a sphere, the formula (2) is subtracted from the formula (1) to obtain the equation for the plane P1 where the two spheres intersect each other, PD1 being positioned on the circle where this plane intersects the sphere represented by the formula (1), this circle having the centre at K₁, the coordinate of which being (a₁, b₁, c₁), and the radius of R₁; and likewise, the formula (4) is subtracted from the formula (3) to obtain the equation for the plane P2 where the two spheres intersect each other, PD2 being positioned on the circle where this plane intersects the sphere represented by the formula (3), this circle having the centre at K₂, the coordinate of which being (a₂, b₂, c₂) and the radius of R₂, the coordinate of the circle centre and the radius of the circle being represented by the formulas:

K _(i)(a _(i) ,b _(i) ,c _(i))=(x _(t1)+(x _(t2) −x _(t1))w _(i) /L, y _(t1)+(y _(t2) −y _(t1))w _(i) /L, z _(t1)+(z _(t2) −z _(t1))w _(i) /L), i=1,2 . . .   (7)

R ₁=√{square root over (d ₁ ² −w ₁ ²)}

R ₂=√{square root over (d ₃ ² −w ₂ ²)}   (8)

where w_(i) represents the distance between LED1 and the plane P_(i) (i=1, 2) where the two spheres intersect each other and L represents the distance between LED1 and LED2.

Second step: as shown in FIG. 3 , two coordinate transformations are performed, where in the first coordinate transformation, an XYZ coordinate system is transformed into an X′Y′Z′ coordinate system, the X′Y′Z′ coordinate system having the circle centre K₁ as its origin of coordinate, the straight line K₁ K₂ as the new Z′ axis, the intersecting line of the circle K₁ and the original XY plane as the new X′ axis; and in the second coordinate transformation, the new X′Y′Z′ coordinate system is transformed into a cylindrical coordinate system, whereupon after the two coordinate transformations, the coordinates of PD1 and PD2 are changed into (R₁, Φ₁, 0) and (R₂, Φ₂, S) respectively, and meanwhile, in the cylindrical coordinate system, the formulas (5) and (6) are represented respectively as:

Φ₁−Φ₂=±arcos M  (9)

β(R ₂ sin Φ₂ −R ₁ sin Φ₁)=−γS  (10)

where M=(R₁ ²+R₂ ²+S²−l²)/(2R₁R₂), the distance between the two circle centres K₁ and K₂ is expressed as S=√{square root over (a²+b²+c²)}, β=−√{square root over ((a²+b²))}/S, γ=c/S, a=a₂−a₁, b=b₂−b₁, c=c₂−c₁.

Third step: Φ₁ and Φ₂ are calculated according to the set of equations (9) and (10), then two inverse coordinate transformations are performed to recover the coordinates of PD1 and PD2 in the XYZ coordinate system through the formulas (11) and (12).

({circumflex over (x)} _(r1) ,ŷ _(r1) ,{circumflex over (z)} _(r1))′=(e _(x) ,e _(y) ,e _(z))(R ₁ cos Φ₁ ,R ₁ sin Φ₁,0)′+(a ₁ ,b ₁ ,c ₁)′  (11)

({circumflex over (x)} _(r2) ,ŷ _(r2) ,{circumflex over (z)} _(r2))′=(e _(x) ,e _(y) ,e _(z))(R ₂ cos Φ₂ ,R ₂ sin Φ₂ ,S)′+(a ₂ ,b ₂ ,c ₂)′  (12)

where

e_(x)=(b/√{square root over (a²+b²)}, −a/√{square root over (a²+b²)}, 0)′, e_(y)=[ac/(S√{square root over ((a²+b²))}), bc/(S√{square root over ((a²+b²))}), −√{square root over ((a²+b²)}/S]′, e_(z)=(a/S, b/S, c/S)′ are the orthogonal basis of the first coordinate transformation.

Fourth step: in the step S3, four sets of solution are obtained by solving Φ₁ and Φ₂, and accordingly, four sets of coordinates of PD1 and PD2 in the XYZ coordinate system are obtained through two inverse coordinate transformations, whereas the actual position includes only one set. As the four sets of solution are spatially symmetrical with respect to the line interconnecting LED1 and LED2, the real solution can be obtained through determination based on the following conditions, including specifically the following steps:

-   -   excluding two sets of solution representing a position above         LED1 and LED2 considering the fact that the position of the         receiver is definitely below LED1 and LED2 in practical         applications; and     -   obtaining the single real solution out of the remaining two sets         of solution by restricting the receiver in movement only on any         side of the plane consisting of LED1, LED2 and the origin.     -   the coordinates of the receiver as (({circumflex over         (x)}_(r1)+{circumflex over (x)}_(r2))/2, (ŷ_(r1)+ŷ_(r2))/2,         ({circumflex over (z)}_(r1)+{circumflex over (z)}_(r2))/2) are         obtained through the single set of solution, in which the         orientation angle η of the receiver has a value expressed by the         formula (13) as:

$\begin{matrix} {\hat{\eta} = {{{sign}\left( {{\hat{y}}_{r2} - {\hat{y}}_{r1}} \right)}{arc}{\cos\left\lbrack \frac{\left( {{\hat{x}}_{r2} - {\hat{x}}_{r1}} \right)}{\sqrt{\left( {{\hat{x}}_{r1} - {\hat{x}}_{r2}} \right)^{2} + \left( {{\hat{y}}_{r1} - {\hat{y}}_{r2}} \right)^{2}}} \right\rbrack}}} & (13) \end{matrix}$

where sign is the sign function.

It is noted that, in the present invention, the range where the terminal is to be positioned is on any side of the plane consisting of LED1, LED2 and the origin.

According to the present invention, accurate 3D positioning and orientation of the terminal can be achieved. Only two LED lamps are used at the transmitting terminal to enable 3D positioning and orientation while providing illumination. A pair of (2 in total) photodetectors is installed on the terminal to receive signals. These two photodetectors are positioned at the same level in the same receiving plane and have a constant and known distance therebetween. According to the present invention, through the TOA (Time of Arrival) principle, the time required for the optical signal to be transmitted from LED1 and LED2 to and received by PD1 and PD2 respectively is measured and the propagation time is multiplied by the speed of light to calculate the distance d₁ between LED1 and PD1, the distance d₂ between LED2 and PD1, the distance d₃ between LED1 and PD2 and the distance d₄ between LED2 and PD2. Positioning is achieved by the calculation steps in the formulas (1) to (13) described above based on the geometrical relationship between the transmitter and the receiver.

The present invention further provides a 3D wireless optical positioning system, including:

LED lamps including a first LED lamp and a second LED lamp with a coordinate of (x_(t1), y_(t1), z_(t1)) and (x_(t2), y_(t2), z_(t2)) respectively, arranged on the ceiling to transmit optical information and provide illumination;

-   -   a receiver provided with a first photodetector and a second         photodetector with a coordinate of ({circumflex over (x)}_(r1),         ŷ_(r1), {circumflex over (z)}_(r1)) and ({circumflex over         (x)}_(r2), ŷ_(r2), {circumflex over (z)}_(r2)) respectively,         arranged on a receiving plane to receive the optical         information;     -   in which the distance between the first photodetector and the         second photodetector is defined as l, meanwhile the first         photodetector and the second photodetector both face upwards on         the receiving plane, the middle point between the first         photodetector and the second photodetector defines the actual         position of the receiving terminal to be predicted, and the         direction from the first photodetector to the second         photodetector defines the orientation of the receiving terminal,         that is, the included angle between the line interconnecting the         first photodetector and the second photodetector and the         positive half of the X axis is the orientation angle η;     -   through the TOA (Time of Arrival) principle, the time required         for the optical signal to be transmitted from the first LED lamp         and the second LED lamp to and received by the first         photodetector and the second photodetector respectively is         measured, and the propagation time is multiplied by the speed of         light to calculate the distance d₁ between the first LED lamp         and the first photodetector, the distance d₂ between the second         LED lamp and the first photodetector, the distance d₃ between         the first LED lamp and the second photodetector and the distance         d₄ between the second LED lamp and the second photodetector; and     -   the actual position and orientation angle of the receiver are         obtained based on the geometrical relationship between the LED         lamps and the receiver in the XYZ coordinate system, the first         photodetector and the second photodetector having a distance         determined as l therebetween and being situated in the same         receiving plane, the range where the receiver is to be         positioned being restricted to be on any side of the plane         consisting of the first LED lamp, the second LED lamp and the         origin and the range of the receiver being restricted to be         below the first LED lamp and the second LED lamp.

EMBODIMENT

To evaluate the performance of the proposed 3D wireless optical positioning method and system, a specific indoor space scenario of a size of 3 m×5 m×3 m (length×width×height) is considered for positioning. The LED lamps are deployed on the ceiling. LED1 has a coordinate of (0, 1.5, 3) and LED2 has a coordinate of (0, 3.5, 3). The distance between the transmitting terminal and the receiving terminal is estimated by using the TOA method. Assuming that a random error Δd is present in distance estimation for d₁, d₂, d₃ and d₄ measured through TOA, and the individual random errors are independent and subject to normal distribution. The mobile terminal is positioned on the receiving plane and has a orientation angle (orientation) that is randomly distributed. Tests are conducted at various positions at an interval of 0.5 m, as shown by the cross line in FIG. 4 . A photodetector PD1 and a photodetector PD2 are installed on the terminal and have a constant distance l therebetween of 0.2 m, and the estimation error Δd has a standard deviation of 0.025 m. As can be seen from the drawings, after the 3D positioning and orientation solution propose by the present invention is adopted, when the receiving plane is at a level of 0.5 m, the 3D positioning error is 9.5 cm and the orientation error is 10.2°; when the receiving plane is at a level of 1 m, the 3D positioning error is 9.1 cm and the orientation error is 7.8°; when the receiving plane is at a level of 1.5 m, the 3D positioning error is 9.0 cm and the orientation error is 10.4°; and three planes at different levels have an overall average 3D positioning error of 9.2 cm and an overall terminal orientation angle error of 9.5°.

FIG. 5 is the diagram of an accumulated distribution function (CDF) of the estimation error for the orientation angle η as the distance 1 between the photodetector PD1 and the photodetector PD2 changes. The level of the receiving plane is at 1 m, the light source 1 is positioned at (0, 1.5, 3) and the light source 2 is positioned at (0, 3.5, 3), and the standard deviation of Δd is 0.025 m. According to the statistical results, when 1 is 0.1 m, 0.2 m, 0.5 m and 1 m, there are 90% orientation angle errors that are less than 4.1°, 2.0°, 0.8° and 0.4° respectively. These results prove the feasibility and reliability of the proposed method.

The embodiments described above are only preferred embodiments listed for fully explaining the present invention, and the scope of protection of the present invention is not limited thereto. Equivalent substitutions or changes made by those skilled in the art on the basis of the present invention shall fall within the scope of protection of the present invention. The scope of protection of the present invention is defined by the claims. 

1. A 3D wireless optical positioning method comprising steps of: arranging LED lamps on a ceiling to transmit optical information and providing illumination and arranging a receiver on a receiving plane to receive the optical information, wherein the LED lamps include a first LED lamp and a second LED lamp with a coordinate of (x_(t1), y_(t1), z_(t1)) and (x_(t2), y_(t2), z_(t2)) respectively, a first photodetector and a second photodetector with a coordinate of ({circumflex over (x)}_(r1), ŷ_(r1), {circumflex over (z)}_(t1)) and ({circumflex over (x)}_(r2), ŷ_(r2), {circumflex over (z)}_(r2)) respectively are arranged on the receiver, a distance between the first photodetector and the second photodetector is defined as l, meanwhile the first photodetector and the second photodetector both face upwards on the receiving plane, a middle point between the first photodetector and the second photodetector defines an actual position of a receiving terminal to be predicted, and a direction from the first photodetector to the second photodetector defines an orientation of the receiving terminal, that is, an included angle between the line interconnecting the first photodetector and the second photodetector and the positive half of the X axis is a orientation angle η; through the Time of Arrival principle, measuring the time required for the optical signal to be transmitted from the first LED lamp and the second LED lamp to and received by the first photodetector and the second photodetector respectively, and multiplying the propagation time by the speed of light to calculate the distance d₁ between the first LED lamp and the first photodetector, the distance d₂ between the second LED lamp and the first photodetector, the distance d₃ between the first LED lamp and the second photodetector and the distance d₄ between the second LED lamp and the second photodetector; and obtaining the actual position and orientation angle of the receiver based on the geometrical relationship between the LED lamps and the receiver in the XYZ coordinate system, the first photodetector and the second photodetector having a distance determined as l therebetween and being situated in the same receiving plane, the receiver being situated below the first LED lamp and the second LED lamp, the range where the receiver is to be positioned being on any side of the plane consisting of the first LED lamp, the second LED lamp and the origin.
 2. The 3D wireless optical positioning method of claim 1, wherein the corresponding photodetector and LED lamp are in time synchronization in calculating the distances d₁-d₄.
 3. The 3D wireless optical positioning method of claim 1, wherein the following equation set is obtained for the distances d₁˜d₄ in the XYZ coordinate system: d ₁ ²=({circumflex over (x)} _(r1) −x _(t1))²+(ŷ _(r1) −y _(t1))²+({circumflex over (z)} _(r1) −z _(t1))²  (1) d ₂ ²=({circumflex over (x)} _(r1) −x _(t2))²+(ŷ _(r1) −y _(t2))²+({circumflex over (z)} _(r1) −z _(t2))²  (2) d ₃ ²=({circumflex over (x)} _(r2) −x _(t1))²+(ŷ _(r2) −y _(t1))²+({circumflex over (z)} _(r2) −z _(t1))²  (3) d ₄ ²=({circumflex over (x)} _(r2) −x _(t2))²+(ŷ _(r2) −y _(t2))²+({circumflex over (z)} _(r2) −z _(t2))²  (4) meanwhile, supplementary equations are obtained as the first photodetector and the second photodetector have a distance determined as l therebetween and are situated in the same receiving plane: l ²=({circumflex over (x)} _(r2) −{circumflex over (x)} _(r1))²+(ŷ _(r2) −ŷ _(r1))²+({circumflex over (z)} _(r2) −{circumflex over (z)} _(r1))²  (5) {circumflex over (z)} _(r2) ={circumflex over (z)} _(r1)  (6) and given the known d₁, d₂, d₃ and d₄ and the formulas (5) and (6) and in combination with the fact that the receiver is situated below the first LED lamp and the second LED lamp and the range where the receiver is to be positioned is on any side of the plane consisting of the first LED lamp, the second LED lamp and the origin, the actual position of the receiver (({circumflex over (x)}_(r1)+{circumflex over (x)}_(r2))/2, (ŷ_(r1)+ŷ_(r2))/2, ({circumflex over (z)}_(r1)+{circumflex over (z)}_(r2))/2) and the orientation angle of the receiver are obtained through solution of: $\begin{matrix} {\hat{\eta} = {{{sign}\left( {{\hat{y}}_{r2} - {\hat{y}}_{r1}} \right)}{arc}{{\cos\left\lbrack \frac{\left( {{\hat{x}}_{r2} - {\hat{x}}_{r1}} \right)}{\sqrt{\left( {{\hat{x}}_{r1} - {\hat{x}}_{r2}} \right)^{2} + \left( {{\hat{y}}_{r1} - {\hat{y}}_{r2}} \right)^{2}}} \right\rbrack}.}}} & (13) \end{matrix}$
 4. The 3D wireless optical positioning method of claim 3, wherein the process of solving the actual position of the receiver (({circumflex over (x)}_(r1)+{circumflex over (x)}_(r2))/2, (ŷ_(r1)+ŷ_(r2))/2, ({circumflex over (z)}_(r1)+{circumflex over (z)}_(r2))/2) and the orientation angle of the receiver specifically comprises steps of: S1: as the formulas (1) and (2) each represent on the physical sense a sphere, subtracting the formula (2) from the formula (1) to obtain an equation for the plane P1 where the two spheres intersect each other, the first photodetector being positioned on a circle where this plane intersects the sphere represented by the formula (1), this circle having a centre at K₁, a coordinate of which being (a₁, b₁, c₁), and a radius of R₁; and likewise, subtracting the formula (4) from the formula (3) to obtain the equation for the plane P2 where the two spheres intersect each other, the second photodetector being positioned on a circle where this plane intersects the sphere represented by the formula (3), this circle having a centre at K₂, a coordinate of which being (a₂, b₂, c₂), and a radius of R₂, a coordinate of the circle centre and a radius of the circle being represented by the formulas: $\begin{matrix} {{{K_{i}\left( {a_{i},b_{i},c_{i}} \right)} = \left( {{x_{t1} + {\left( {x_{t2} - x_{t1}} \right)\frac{w_{i}}{L}}},{y_{t1} + {\left( {y_{t2} - y_{t1}} \right)\frac{w_{i}}{L}}},{z_{t1} + {\left( {z_{t2} - z_{t1}} \right)\frac{w_{i}}{L}}}} \right)},{i = 1},{2\ldots}} & (7) \end{matrix}$ $\begin{matrix} \begin{matrix} {R_{1} = \sqrt{d_{1}^{2} - w_{1}^{2}}} \\ {R_{2} = \sqrt{d_{3}^{2} - w_{2}^{2}}} \end{matrix} & (8) \end{matrix}$ where w_(i) represents a distance between the first LED lamp and a plane P_(i) (i=1, 2) where the two spheres intersect each other and L represents a distance between the first LED and the second LED lamp; S2: performing two coordinate system transformations, where in a first coordinate transformation, an XYZ coordinate system is transformed into an X′Y′Z′ coordinate system, the X′Y′Z′ coordinate system having a circle centre K₁ as its origin of coordinate, a straight line K₁ K₂ as a new Z′ axis, an intersecting line of the circle K₁ and the original XY plane as a new X′; and in the second coordinate transformation, a new X′Y′Z′ coordinate system is transformed into a cylindrical coordinate system, whereupon after the two coordinate transformations, the coordinates of the first photodetector and the second photodetector are changed into (R₁, Φ₁, 0) and (R₂, Φ₂, S) respectively, and meanwhile, in a cylindrical coordinate system, the formula (5) and the formula (6) are represented respectively as: Φ₁−Φ₂=±arcos M  (9) β(R ₂ sin Φ₂ −R ₁ sin Φ₁)=−γS  (10) where M=(R₁ ²+R₂ ²+S²−l²)/(2R₁R₂), a distance between the two circle centres K₁ and K₂ is expressed as S=√{square root over (a²+b²+c²)}, β=−√{square root over ((a²+b²))}/S, γ=c/S, a=a₂−a₁, b=b₂−b₁, c=c₂−c₁; and S3: calculating Φ₁ and Φ₂ according to the formula (9) and the formula (10), then performing two inverse coordinate transformations to recover the coordinates of the first photodetector and the second photodetector in the XYZ coordinate system through the formulas (11) and (12): ({circumflex over (x)} _(r1) ,ŷ _(r1) ,{circumflex over (z)} _(r1))′=(e _(x) ,e _(y) ,e _(z))(R ₁ cos Φ₁ ,R ₁ sin Φ₁,0)′+(a ₁ ,b ₁ ,c ₁)′  (11) ({circumflex over (x)} _(r2) ,ŷ _(r2) ,{circumflex over (z)} _(r2))′=(e _(x) ,e _(y) ,e _(z))(R ₂ cos Φ₂ ,R ₂ sin Φ₂ ,S)′+(a ₂ ,b ₂ ,c ₂)′  (12) where e_(x)=(b/√{square root over (a²+b²)}, −a/√{square root over (a²+b²)}, 0)′, e_(y)=[ac/(S√{square root over ((a²+b²))}), bc/(S√{square root over ((a²+b²))}), −√{square root over ((a²+b²)}/S]′, e_(z)=(a/S, b/S, c/S)′ are the orthogonal basis of the first coordinate transformation.
 5. The 3D wireless optical positioning method of claim 4, wherein in the step S3, four sets of solution are obtained by solving Φ₁ and Φ₂, and accordingly, four sets of coordinates of the first photodetector and the second photodetector in the XYZ coordinate system are obtained through two inverse coordinate transformations, whereas the actual position includes only one set, and as the four sets of solution are spatially symmetrical with respect to the line interconnecting the first LED lamp and the second LED lamp, the real solution is obtained through determination based on the following conditions, including steps of: S4: excluding two sets of solution representing a position above the first LED lamp and the second LED lamp considering the fact that the position of the receiver is below the first LED lamp and the second LED lamp; S5: obtaining a single real solution out of the remaining two sets of solution by restricting the receiver in movement on any side of the plane consisting of the first LED lamp, the second LED lamp and the origin; S6: obtaining the coordinate of the receiver as (({circumflex over (x)}_(r1)+{circumflex over (x)}_(r2))/2, (ŷ_(r1)+ŷ_(r2))/2, ({circumflex over (z)}_(r1)+{circumflex over (z)}_(r2))/2) through the single set of solution, in which the orientation angle η of the receiver has a value expressed by the formula (13) as: $\begin{matrix} {\hat{\eta} = {{{sign}\left( {{\hat{y}}_{r2} - {\hat{y}}_{r1}} \right)}{arc}{\cos\left\lbrack \frac{\left( {{\hat{x}}_{r2} - {\hat{x}}_{r1}} \right)}{\sqrt{\left( {{\hat{x}}_{r1} - {\hat{x}}_{r2}} \right)^{2} + \left( {{\hat{y}}_{r1} - {\hat{y}}_{r2}} \right)^{2}}} \right\rbrack}}} & (13) \end{matrix}$ where sign is the sign function.
 6. The 3D wireless optical positioning method of claim 1, wherein the number of the LED lamps is defined depending on the region where they are to be positioned.
 7. A 3D wireless optical positioning system comprising: LED lamps including a first LED lamp and a second LED lamp with a coordinate of (x_(t1), y_(t1), z_(t1)) and (x_(t2), y_(t2), z_(r2)) respectively, arranged on the ceiling to transmit optical information and provide illumination; a receiver provided with a first photodetector and a second photodetector with a coordinate of ({circumflex over (x)}_(r1), ŷ_(r1), {circumflex over (z)}_(r1)) and ({circumflex over (x)}_(r2), ŷ_(r2), {circumflex over (z)}_(r2)) respectively, arranged on a receiving plane to receive the optical information; wherein a distance between the first photodetector and the second photodetector is defined as l, meanwhile the first photodetector and the second photodetector both face upwards on the receiving plane, a middle point between the first photodetector and the second photodetector defines an actual position of a receiving terminal to be predicted, and a direction from the first photodetector to the second photodetector defines an orientation of the receiving terminal, that is, an included angle between a line interconnecting the first photodetector and the second photodetector and the positive half of the X axis is a orientation angle η; through the Time of Arrival principle, the time required for the optical signal to be transmitted from the first LED lamp and the second LED lamp to and received by the first photodetector and the second photodetector respectively is measured, and the propagation time is multiplied by the speed of light to calculate a distance d₁ between the first LED lamp and the first photodetector, a distance d₂ between the second LED lamp and the first photodetector, a distance d₃ between the first LED lamp and the second photodetector and a distance d₄ between the second LED lamp and the second photodetector; and an actual position and orientation angle of the receiver are obtained based on a geometrical relationship between the LED lamps and the receiver in the XYZ coordinate system, the first photodetector and the second photodetector having a distance determined as l therebetween and being situated in the same receiving plane, the receiver being situated below the first LED lamp and the second LED lamp, the range where the receiver is to be positioned being on any side of the plane consisting of the first LED lamp, the second LED lamp and the origin.
 8. The 3D wireless optical positioning system of claim 7, wherein the first LED lamp and the second LED lamp and the first photodetector and the second photodetector have synchronized operation time.
 9. The 3D wireless optical positioning system of claim 7, wherein a set of equations for the distances d₁ to d₄ in the XYZ coordinate system is obtained based on the geometrical relationship between the LED lamps and the receiver in the XYZ coordinate system: d ₁ ²=({circumflex over (x)} _(r1) −x _(t1))²+(ŷ _(r1) −y _(t1))²+({circumflex over (z)} _(r1) −z _(t1))²  (1) d ₂ ²=({circumflex over (x)} _(r1) −x _(t2))²+(ŷ _(r1) −y _(t2))²+({circumflex over (z)} _(r1) −z _(t2))²  (2) d ₃ ²=({circumflex over (x)} _(r2) −x _(t1))²+(ŷ _(r2) −y _(t1))²+({circumflex over (z)} _(r2) −z _(t1))²  (3) d ₄ ²=({circumflex over (x)} _(r2) −x _(t2))²+(ŷ _(r2) −y _(t2))²+({circumflex over (z)} _(r2) −z _(t2))²  (4), meanwhile, supplementary equations are obtained as the first photodetector and the second photodetector have a distance determined as l therebetween and are situated in the same receiving plane: l ²=({circumflex over (x)} _(r2) −{circumflex over (x)} _(r1))²+(ŷ _(r2) −ŷ _(r1))²+({circumflex over (z)} _(r2) −{circumflex over (z)} _(r1))²  (5) {circumflex over (z)} _(r2) ={circumflex over (z)} _(r1)  (6); and given the known d₁, d₂, d₃ and d₄ and the formulas (5) and (6) and in combination with the fact that the receiver is situated below the first LED lamp and the second LED lamp and the range where the receiver is to be positioned is on any side of the plane consisting of the first LED lamp, the second LED lamp and the origin, the actual position of the receiver (({circumflex over (x)}_(r1)+{circumflex over (x)}_(r2))/2, (ŷ_(r1)+ŷ_(r2))/2, ({circumflex over (z)}_(r1)+{circumflex over (z)}_(r2))/2) and the orientation angle of the receiver are obtained through solution of: $\begin{matrix} {\hat{\eta} = {{{sign}\left( {{\hat{y}}_{r2} - {\hat{y}}_{r1}} \right)}{arc}{{\cos\left\lbrack \frac{\left( {{\hat{x}}_{r2} - {\hat{x}}_{r1}} \right)}{\sqrt{\left( {{\hat{x}}_{r1} - {\hat{x}}_{r2}} \right)^{2} + \left( {{\hat{y}}_{r1} - {\hat{y}}_{r2}} \right)^{2}}} \right\rbrack}.}}} & (13) \end{matrix}$
 10. The 3D wireless optical positioning system of claim 9, wherein the first LED lamp, the second LED lamp, the first photodetector and the second photodetector are provided with a time synchronization device. 